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ABSTRACT 

We report the detection of surprisingly strong HCN(J=6^5), HNC(J=6^5), and HCO+(J=6^5) 
emission in the host galaxy of the z=3.91 quasar APM 08279+5255 through observations with the 
Combined Array for Research in Millimeter- wave Astronomy (CARMA). HCN, HNC, and HCO + are 
typically used as star formation indicators, tracing dense molecular hydrogen gas [n(H2) > 10 5 cm -3 ] 
within star- forming molecular clouds. However, the strength of their respective line emission in the 
J=6^5 transitions in APM 08279+5255 is extremely high, suggesting that they are excited by another 
mechanism besides collisions in the dense molecular gas phase alone. We derive J=6— ?>5 line lumi- 
nosities of L^ CN = (4.9 ± 0.6), L^ NC = (2.4 ± 0.7), and L' HCO+ = (3.0 ± 0.6) x 10 10 /i^ 1 K kms" 1 pc 2 
(where is the lensing magnification factor), corresponding to V ratios of ~0. 23-0. 46 relative to 
CO(J=l— >-0). Such high line ratios would be unusual even in the respective ground-state (J=l— »>0) 
transitions, and indicate exceptional, collisionally and radiatively driven excitation conditions in the 
dense, star-forming molecular gas in APM 08279+5255. Through an expansion of our previous mod- 
eling of the HCN line excitation in this source, we show that the high rotational line fluxes are caused 
by substantial infrared pumping at moderate opacities in a ~220K warm gas and dust component. 
This implies that standard M^ ense /L f conversion factors would substantially overpredict the dense 
molecular gas mass Md en se- We also find a HCN(J=6^5)/HCN(J=5^4) L' ratio greater than 1 
(1.36+0.31) - however, our models show that the excitation is likely not 'super-thermal', but that the 
high line ratio is due to a rising optical depth between both transitions. These findings are consistent 
with the picture that the bulk of the gas and dust in this source is situated in a compact, nuclear 
starburst, where both the highly active galactic nucleus and star formation contribute to the heating. 
Subject headings: galaxies: active — galaxies: starburst — galaxies: formation — galaxies: high- 
redshift — cosmology: observations — radio lines: galaxies 



1. INTRODUCTION 

Over the past decade, great progress has been made in 
understanding the conditions for star formation in gas- 
rich galaxies out to the highest redshifts. Molecular gas, 
the prospective fuel for star formation, is now detected in 
more than 70 galaxies at z > 1, allowing to compare the 
properties of star-forming environments among different 
galaxy popul ations in the early universe (see Solomon & 
Vanden Bout 2005 for a review). These detections are al- 
most exclusively being obtained in rotational transitions 
of CO, which (due to its relatively low critical density 
of ^crit(H2) — 300 cm -3 ) is a good proxy for the total 
amount of molecular gas in a galaxy. 

More focused studies of the dense molecular gas found 
in the star-forming cores of molecular clouds typically 
employ observations of high dipole moment, high critical 
density (n cr i t (H2) > 10 4 cm -3 ) molecules such as HCN, 
HCO + , and HNC, both in nearby galaxies and out to 
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high z (e.g., Gao & Solomon 2004; Rieche rs et al. I2006al 
l2007at_Gao et al. 120071 : Baan et al. 120081 : Gracia-Carpio 
et al. 120081 ) . 

We here aim to study, for the first time, the dense, 
star-forming molecular gas excitation in a high-z galaxy. 
Constraints on the excitation of gas at very high den- 
sities and the physical mechanisms (collisions vs. other 
channels) responsible are crucial to understand in more 
detail how excitation may influence scaling relations be- 
tween the dense gas content and star formation rate of 
galaxies back to early cosmic times. To disentangle exci- 
tation effects from other phenomena (such as the chem- 
ical composition of the gas), it is crucial to target mul- 
tiple diagnostic lines in a well-studied, key system. The 
target of this study is the extremely luminous z=3.91 
quasar APM 08279+5255, which has been studied com- 
prehensively in CO and o ther d iagnostics (see, e.g., WeiB 
et al. 2007; Riechers et al. 120091 for details). In particular, 
it is one out of only two high-z galaxies in which mul- 
tiple dense molecular gas tracers were detected to date 
(the other being the Clover leaf quasa r at z=2.56; Bar- 



2005 



2006; 



vamis et al. 1997; Solomon et al. 2003; Wagg et al. 
Riechers et al. I2006al l2007ri Garcia-Burillo et al 
Guelin et al. 120071 1. 

In this paper, we report the detection of 
HCN(J=6^5), HCO+(J=6^5), and HNC(J=6^5) 
emission in the quasar host galaxy of APM 08279+5255 
(2=3.91), using the Combined Array for Research in 
Millimeter- wave Astronomy (CARMA). These observa- 



2 



Riechers et al. 



109.0 



108.5 



108.0 



"obs 



111.0 



[GHz] 
110.5 



106.0 105.5 105.0 104.5 



8 R 



E 



CO 

c 

"D 

X 
3 



6 - 



4 - 



HCN * & HC0 + ( J=6-5) 



1 1 1 I 1 1 1 1 I 1 1 1 
LSB continuum 




-1000 



1000 -1000 1000 

Velocity offset [km/s] 



-2000 



2000 



Fig. 1. CARMA spectra of the HCN(J=6^5), HCO+(J=6^5) (left), HNC(J=6^5) (middle), and continuum emission at 2.76, 2.71, 
and 2.84 mm (left to right) toward APM 08279+5255 at a resolution of ~85 kms -1 (31.25 MHz). The velocity scales are relative to the 
tuning frequencies of 108.611GHz (left) and 110.751GHz (middle) and the central frequency of both LSBs at 105.391GHz (right). The 
rms per velocity bin are 0.95, 0.88, and ~0.65-0.92 mJy (left to right). The solid curves show simultaneous Gaussian fits to the line and 
continuum emission where applicable. 



tions represent the first extragalactic detections of such 
high- J lines of the dense gas tracers HCN, HCO + , and 
HNC, and significantly constrain the physical properties 
of the dense gas in the star-forming regions of this 
distant galaxy. We use a concordance, flat ACDM 
cosmology throughout, with Hq= 71 kms" 1 M pc -1 , 
^ M =0.27, and ft A =0.73 (Spergel et al. [2QQ3L l2007h . 

2. OBSERVATIONS 

We used CARMA to observe the HCN( J=6^5) rest 
= 531.7164GHz), HCO+(J=6^5) (535.0618 GHz), 
and HNC(J=6^5) (543.8974 GHz) transition lines 
(-560 /im) toward APM 08279+5255. At ^=3.911, these 
lines are redshifted to 108.270, 108.952, and 110.751 GHz 
(~ 2.8mm). The target was observed with 14 or 15 an- 
tennas (corresponding to 91 or 105 baselines per an- 
tenna configuration) for 14 tracks in March and June 
2008 (setup 1; HCN/HCO+, observed simultaneously) 
and 23 tracks between February and July 2009 (setup 2; 
HNC), amounting to a total observing time of 148 hr. 
This results in 37 hr on source time for HCN/HCO + and 
52 hr for HNC. All HCN/HCO+ observations were car- 
ried out in D array (11-148 m baselines), while HNC ob- 
servations were carried out in the C, D and E arrays 
(6-185 m baselines after flagging). Observations before 
June 2008 were carried out with the previous genera- 
tion 3 mm receivers, and observations after March 2008 
were carried out with the new generation 3 mm receivers 
(which offer improved noise temperatures, tuning range, 
and stability). 

Weather conditions scaled between acceptable and 
excellent for observations at 3 mm wavelengths. For 
the HCN/HCO + observations, typical median phase 
rms values were 275 /am (normalized to a 100 m base- 
line, measured at 45° elevation), median optical depths 
were T23ogh z =0.24, and median precipitable water va- 
por columns were 3.9 mm. For the HNC observations, 
typical median phase rms values were 300 (C arrayjH, 



235 (D array), and 395 /im (E array), median optical 
depths were T23ogh z =0.38, and median precipitable wa- 
ter vapor columns were 5.5 mm. The nearby source 
J0818+423 (distance to APM 08279+5255: 10.6°) was 
observed every 20 minutes for secondary amplitude and 
phase calibration. The strong calibrator sources J0423- 
013, J0927+390, 3C111, 3C84, and 3C273 were ob- 
served at least once per track for bandpass and secondary 
flux calibration. Absolute fluxes were bootstrapped rel- 
ative to Mars, Uranus, or 3C 84 (when no planet was 
available). Pointing was performed at least every 2-4 hr 
on nearby sources, using both radio and optical modes. 
The resulting total calibration is estimated to be accu- 
rate within —15% (see Appendix). 

The 3 mm receivers were tuned between the HCN and 
HCO+ lines at 108.611GHz (setup 1) and at the red- 
shifted HNC frequency of 110.751 GHz (setup 2), both in 
the upper sideband (USB). The intermediate frequencies 
were chosen to be 1.79 GHz for the HCN/HCO+ obser- 
vations (optimized for the previous generation receivers) 
and 2.50 GHz for the HNC observations (optimized for 
the new generation receivers), centering the lower side- 
bands (LSBs) at 105.031 and 105.751GHz. Three bands 
with 15 channels of 31.25 MHz (~85 kms -1 ) width each 
were centered on the tuning frequencies. The bands were 
overlapped by 2 channels to improve calibration of the 
correlated dataset, leading to an effective bandwidth of 
1281.25MHz (-3500 kms -1 ) per sideband. This com- 
fortably covers the HCN, HCO + , and HNC emission 
lines as well as the 2.8 mm (rest-frame 560 jam) con- 
tinuum of the source. It also provides coverage of 
the redshifted HC 3 N J=57^56, 59^58, and 60^59 
lines (^ rest = 518.1897, 536.3420, and 545.4170 GHz; 
^obs = 105.516, 109.212, and 111.060GHz), as well as 
the SiO(J=12->ll) line (^ rest = 520.8782 GHz; is ohs = 
106.064 GHz). In particular, HC 3 N( J=57^56) lies in a 
spectral range that is covered by both LSB frequency 
setups. 



7 Only one 3 hr track was observed in C array, and all baselines 



>185m were discarded due to the relatively high phase noise. 
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TABLE 1 

Measured line fluxes and luminosities in APM 08279+5255. 





I 


V 




[Jy km s ~ 1 ] 


[10 10 /i -1 K kms-ipc 2 ] 


HCO+(J=6^5) 


1.01 ± 0.19 


3.0 ± 0.6 


HCN(J=6^5) 


1.65 ± 0.19 


4.9 ± 0.6 


HNC(J=6^5) 


0.86 ± 0.24 


2.4 ± 0.7 


HC 3 N(J=57^56) a 


<0.34 


<1.0 


HC 3 N(J=59^58) 


<0.52 


<1.5 


HC 3 N(J=60^59) 


<0.49 


<1.4 


SiO(J=12->41) 


<0.44 


<1.3 



Note. — All limits are 3a, extracted over a linewidth of 
400 kms -1 . Luminosities are derived as described by Solomon etal. 
(fl992h : L'fKkms^pc 2 ] = 3.25 X 10 7 X J X u~^ s x X (1 + z)~ 3 , 
where / is the velocity integrated line flux in Jy kms -1 , Dl is the 
luminosity distance in Mpc (z = 3.911), and u bs is the observed fre- 
quency in GHz. The given luminosities are not correc ted for the lensing 
magnification factor of /2l=4.2 (Riechers et al. |2009f) . 

a Covered by both frequency settings. 

For data reduction and analysis, the MIRIAD package 
was used. The final plots were created with the GILDAS 
package. All data were imaged using 'natural' weight- 
ing, resulting in synthesized clean beams of 4.5" x 3.6" 
(HCN/HCO+), and 5.3"x4.6" (HNC; -29 and 35kpc at 
z=3. 91). T he CO size of the source is <0.4" (Riechers 
et al. 120091 ), i.e., it is unresolved at the resolution of our 
observations. The final rms values for the HCN/HCO + 
and HNC observations are 0.15 and 0.13 mJy beam -1 
over the full USB bandpass (1281.25 MHz), and 0.95 and 
0.88 mJy beam -1 per 31.25 MHz (—85 kms -1 ) channel. 
The rms noise in the LSB is by 3% (HCN/HCO+) and 
15% (HNC) better compared to the USB due to atmo- 
spheric transmission. 

3. RESULTS 

3.1. Emission Lines of High- Density Gas Tracers 

We have detected emission from the J=6^5 HCN, 
HCO + and HNC emission lines toward the z=3.91 quasar 
APM 08279+5255. The left panel of Figure ffl shows the 
HCN and HCO+ spectrum at 86 kms -1 (31.25 MHz) 
resolution. Both lines are detected simultaneously on 
top of strong 2.76 mm continuum emission. The middle 
panel of Figure [1] shows the HNC spectrum at similar 
(85 kms -1 ; 31.25 MHz) resolution. The line is also de- 
tected on top of strong continuum emission at similar 
wavelength (2.71mm). Figure [2] shows the HCN and 
HCO + velocity channel maps at the same velocity res- 
olution and coverage as the spectrum in Fig. [TJ At an 
rms of 0.95 mJy beam -1 , APM 08279+5255 is detected 
at 8<j and 6a significance in the peak channels of the 
HCN and HCO + lines. Both emission lines are detected 
over 10 channels on top of the continuum, showing a de- 
cline in strength toward the linewings, as expected. Also, 
marginal excess flux is detected at the peak position of 
the redshifted HCaN( J=59— ^58) line, which however is 
not formally detected. Figure [3] shows the HNC veloc- 
ity channel maps at the same velocity resolution as the 
spectrum in Fig. [T] (central 12 channels). At an rms of 
0.88 mJy beam -1 , the source is detected at 6a signifi- 
cance in the peak channel of the HNC line, also showing 
a decline in flux towards the line wings. 

From simultaneous Gaussian fitting to the profiles 
of the HCN(J=6^5) and HCO+(J=6^5) lines and 
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TABLE 2 

Measured continuum fluxes in APM 08279+5255. 



^obs 

[GHz] 


^obs 

[mm] 


^rest 

H 


[mJy] 


110.7 


2.71 


551 


2.17+0.19 


108.6 


2.76 


562 


2.08+0.22 


105.4 a 


2.84 


579 


2.04+0.08 


107.5 b 


2.79 


568 


2.08+0.07 



a Combined LSB data from both setups. 
b Averaged over all continuum data. 

the underlying continuum emission, we derive HCN 
and HCO + line peak flux densities of 4.04+0.61 and 
2.48+0.56 mJy and FWHM velocities of 385+36 and 
385+60 kms -1 , respectively (in good agreement with 
the 400+ 40 kms" 1 FWHM of HCN J=5^4; WeiB et 
al. [2007). This corresponds to integrated HCN and 
HCO+ line fluxes of 1.65+0.19 and 1.01+0.19 Jy kms" 1 . 
A simultaneous fit to the HNC(J=6^5) line and un- 
derlying continuum emission yields a HNC line peak 
flux density of 1.57+0.51 mJy and a line FWHM of 
520+160 kms -1 (also consistent with the FWHM of 
HCN J=5— ?>4 within the errors). This corresponds to an 
integrated HNC line flux of 0.86+0.24 Jy kms -1 . From 
the line peak velocities, we derive formal line redshifts 
of ^=3.9126+0.0011 for HCN/HCO+ (simultaneous fit), 
and ^=3.9144+0.0011 for HNC, which are consistent 
with those derived from other molecular emission lines 
withi n the errors (e.g., Riechers et al. 2006b; WeiB et al. 
l2007h . 

We also searched for HC 3 N J=57^56, 59^58, and 
60— >-59, as well as SiO(J=12— >T1) emission within the 
covered spectral range, including the LSBs (combined 
spectrum shown in the right panel of Fig. [TJ. Marginal 
excess flux is seen close to the peak position of the red- 
shifted HCaN(J=57— )>56) line in the overlap region of 
both LSB setups, which we however consider not de- 
tected in the following. We place 3a limits of 0.34, 
0.52, 0.49, and 0.44 Jy kms -1 on the integrated fluxes 
from these lines, extracted over a fixed velocity range of 
400 kms -1 (see Table [TJ. We also attempted to stack all 
three HC3N lines, which results in no clear signal above 
the formal 3a limit of 0.26 Jy kms -1 . We however note 
that due to the uncertainties in extraction, this 'stacked' 
limit has to be treated with caution. 

3.2. Millimeter Continuum Emission 

As mentioned above, we have detected ~2.8mm (rest- 
frame ~560 /im) continuum emission toward the host 
galaxy of the 2=3.91 quasar APM 08279+5255. Emission 
was detected at high signal-to-noise in each sideband of 
the two frequency setups. Thus, continuum fluxes were 
extracted by fitting a two-dimensional, elliptical Gaus- 
sian to the source in the u — v plane for each sideband and 
frequency setting, excluding ranges where line emission 
was detected. The individual values are listed in Table [2j 
and are fully consistent with the continuum fluxes ob- 
tained from the simultaneous line/continuum fits to the 
spectra as outlined above. A combination of all measure- 
ments yields an average continuum flux of 2.08+0.07 mJy 
at 2.79 mm, consistent with the spectral energy distribu- 
tion (SED) of the source (Riechers et al. 120091 ). 
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Fig. 2. — Channel maps of the data shown as a spectrum in Fig.[I](left) at the same velocity resolution. Frequencies increase with channel 
number from left to right, i.e., red to blue. The peak channels of the redshifted HCN, HCO + , and HC3N lines are indicated. Contours are 
shown at (-3, -2, 2, 3, 4, 5, 6, 7, 8) a (1 a = 0.95 mjy beam -1 ). The beam size (4.5" X3.6") is shown in the bottom left corner. The cross 
indicates a position of 08 h 31 m 41 s .70, +52°45 / 17 // .5 between the CO lens images (Riechers et al. [2009h . 



3.3. Line Luminosities and Ratios 
From the line intensities, we derive line luminosities 
and limits of ^hcn(6-5) = ( 4 - 9 ± °- 6 )> L hco+(6-5) = 
(3.0 ±0.6), L' wc(6 _ 6) = (2.4 ±0.7), ^ C3 N ( 57-56) < 

1.0, £ H C 3 N(59-58) < ^HC 3 N(60-59) < anC * 

L SiO(i2-ii) < L3 x lO^/iL 1 ^ kms^pc 2 (Tab. dj not 
corrected for t he len sing magnification factor of /iL=4.2; 
Riechers et al. [2009T). 



This corresponds to J=6^5 V ratios of HCN/CO = 
0.46±0.07, HCO+/CO = 0.28±0.06, and HNC/ CO = 
0.23±0.07 (relative to CO J=1^0; Riechers et al. 120091 ) . 
Such high line ratios would be unusual even in the re- 
spective ground-state ( J=l— K)) transitions, and indicate 
exceptional dense gas excitation conditions in this source. 

We also find a HCN/HCO+ J=6^5 V ratio of 
1.7±0.4, which is higher than the HCN/HCO+ J=5^4 
V ra tio of 1.0±0.2 (WeiB et al. 120071 : Garcia-Burillo et al. 
I2006f ). Such a steeply rising HCN/HCO+ line ratio with 



HCN, HNC, and HCO+(J=6^5) in APM 08279+5255 0=3.91) 



4- 



O" 



+ 



52°45'30" 



52 45'15" 



52 o 45'00" 



=^6-5)' 



J . /O . L 



_l i i a r 



Fig. 3. — Channel maps of the central 12 channels of the data shown as a spectrum in Fig. ^ (middle) at the same velocity resolution. 
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Contours are shown at (-3, -2, 2, 3, 4, 5, 6) a (la = 0.88 mJy beam -1 ). The beam size (5. 3" x4.6 // ) is shown in the bottom left corner. 
The cross indicates the same position as in Fig. [2] 



J is not expected, as the densities and temperatures to 
collisionally excite HCN and HCO + are similar. We note 
that the J=5^4 lines were not observed simultaneously, 
so the error bars of their ratio may be underestimated 
due to the uncertainties in their relative flux scales. In 
addition, the difference in line ratio with J may indicate 
that HCaN( J=49^48) line emission contributes signifi- 
cantly to the HCO + (J=5^4) flux reported by Garcia- 
Burillo et al., as these lines are blended. The HC3N limits 
derived above would be consistent with up to a 30-40% 
contamination (extrapolated from these higher- J transi- 
tions) of the HCO + (J=5^4) flux. Note that such high 
HCsN/HCO + line ratios are actually observed in nearby 
infrared-lu minous galaxies (albeit in lower- J lines; e.g., 
Aalto et al. l2007at ). 

Furthermore, we find HNC(J=6^5)/HCN(J=6^5) 
= 0.50±0.15 and HNC(J=6^5)/HCO+(J=6^5) = 
0.82±0.27. The HNC(J=5^4) line was observed to- 
ward APM 08279+5255, bu t it is strongly blended with 
CN(7V=4^3) (Guelin et al. l2007h . The single Gaussian 
fit by Guelin et al. to the blend of these lines can be trans- 
lated to upper limits of HNC( J=5^4)/HCN( J=5^4) < 
1.12 and HNC(J=5^4)/HCO+(J=5^4) < 1.15. Guelin 
et al. estimate that the CN line contributes ~l/3 to 
the total flux they measure. Assuming their decompo- 
sition is correct, we find HNC( J=5^4)/HCN( J=5^4) 
= 0.64±0.32 and HNC(J=5-»4)/HCO+(J=5->4) = 
0.66+0.33. Given the similar excitation densities and 
temperatures of HNC compared to HCN and HCO + , we 
assume that the J=5^4 and 6— ?>5 line ratios should be 
comparable. This assumption is consistent with the ob- 
served HNC/HCN ratios if the CN contribution to the 
HNC( J=5^4) V is as suggested by Guelin et al. It also 
is consistent with the observed HNC/HCO + ratios, in 
particular if the HC3N contribution to HCO + ( J=5— ?>4) 
is comparable to the CN contribution to HNC( J=5— ?>4). 
However, given the remaining uncertainties, observations 
of another CN transition would be desirable to test this 
scenario (see also discussion below). 

Interestingly, we find a HCN( J=6->5)/HCN( J=5->4) 
V ratio of 1.36+0.31. A ratio of greater than 1 is not ex- 
pected for optically thick, collisionally excited emission 
that is thermalized or subthermal. However, such high 
L' ratios have also been found for HNC in nearby (ultra-) 



luminous infrared galaxies ((U)LIRGs ), suc h as Arp220 
(albeit in lower- J lines; Aalto et al. 2007b). Moreover, 
extreme gas densities of n(H2) > 10 6 cm -3 would be re- 
quired to highly excite the J=6^5 transitions of HCN, 
HNC and HCO+ collisionally (n crit (H 2 ) > 10 8 cm- 3 ). 

4. HCN LINE EXCITATION MODELING 

To understand the unusual dense gas excitation in 
APM 08279+5255 in more detail, we have carried out 
Large Velocity Gradient (L VG) m odels of the HCN line 
excitation (see WeiB et al. 120071 for initial study). Our 
models take both collisional and radiative excitation via 
infrared (IR) pumping into acc ount, using the HCN col- 
lision rates from Schoier et al. (|2005f ). and including the 
first ro- vibrational bending m ode at 14.0 /mi for pumping 
(i/ 2 = l; e.g., Thorwirth et al. l2003f ). We adopted a HCN 
abundance per velocity gradient of [H CN]/( d^/dr) = 
1 x 10_ 9 pc(kms- 1 )- 1 (Heifer & Blitz [l997|; Wang et 
al. 120041 ). All parameters are required to also repro duce 
the CO excitation of the source (WeiB et al. 2007), as- 
suming [HCN/CO]=10 -4 , and are required to be con- 
sistent with the dust spectr al ene rgy distribution (WeiB 
et al. 120071 : Riechers et al. I2009f ). The best solutions 
were obtained for a spherical, single-component model 
with kinetic temperatures of Tki n =220K, gas densities of 
n gas =10 4 - 2 cm~ " 
of Ti R (=T dust ) 

of IRff =0.3-0.7 (Fig. IU where solutions with higher IRff 
prefer slightly lower T^ n to stay consistent with the ob- 
served size of the source in CO emission; Riechers et al. 
l2009h . 

Due to the high intensity of the IR-pumping field (as 
indicated by IRff) in these models, HCN(J=5— ?>4) and 
HCN( J=6^5) are only moderately optically thick (typ- 
ical optical depths of T5_4=0.25 and T6-5=0.4), and the 
optical depth still rises with J between these lines. This 
implies that the excitation of HCN is not 'super-thermal' 
(as in a 'true' population inversion between energy lev- 
els), but that the high HCN L' ratio is merely an optical 
depth effect. The HCN excitation is dominated by IR 
pumping for all models that fit the data. Interestingly, 

8 IRff describes the solid angle fraction of the gas that is exposed 
to the IR field, which in turn is described by a grey-body spectrum. 
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Fig. 4. — HCN excitation ladder (spectral line energy distri- 
bution; points) and LVG models (lines) for APM 08279+5255, 
accounting for both collisional and radiative excitation. The 
HCN(J=5^4) data point is from WeiB et al. ((20071) . The models 
give n g as=10 4 - 2 cm- 3 , T kin =T IR =220 K, and infrared field filling 
factors of IRff=0.3-0.7 (shown in steps of 0.1). 

the models imply that the HCN emission cannot dom- 
inantly arise from a relativel y cold , dense gas and dust 
component (see WeiB et al. 2007 for multi-component 
models). This means that standard conversion factors 
c^hcn = Mdense/^HCN fr° m HCN lin e lum inosity to 
dense gas mass (e.g., Gao & Solomon l2004f ) will sub- 
stantially overpredict the dense molecular gas mass in 
this system. Also, the models do not require unusually 
high relative HCN abundances, contrary to previous sug- 
gestions based on a narrower exploration of the param- 
eter space in excitation conditions (Garcia-Burillo et al. 
l2006h . 

5. DISCUSSION 

5.1. HCO+ and HNC Excitation 

Due to the remaining uncertainties in the 
HCO+(J=5^4) and HNC(J=5->4) line intensities 
(see discussion above), it is currently not possible 
to properly constrain similar models for the HCO + 
and HNC line excitation. However, due to similar 
ro- vibrational bending modes {y^ — 1 at 12.1 and 
21.6 /im), it is likely that similar conclusions hold for 
the HCO + and HNC excitation. The bending modes 
of HCN, HCO+ and HNC all lie close to the peak of 
the dust SEP of A PM 08279+5255 (WeiB et al. 120071 
Riechers et al. l2009f ). which would make such a scenario 
plausible. Also, the organic compound HC3N has 
several bending mode s in this wavelength regime (e.g., 
Wyrowski et al. 1999), which would make a substantial 
HC 3 N(J=49^48) contribution to HCO+(J=5^4) 
plausible. Luminous, high-J HC3N emission consistent 
with radiative excitation scenarios is o bserved in nearby 
galaxies like NGC4418 (Aalto et al. I2007al ). On the 
other hand, the lowest ro-vibrational transition of the 
diatomic molecule CN lies at 4.9 /im (see also Guelin et 
al. 120071 ). which lies substantially beyond the peak of 



the SED. No strong 4.9 /im (observed-frame 24 /im) CN 
absorption is detected towards the bright mid-infrared 
continuum of APM 08279+5255 in deep Spitzer Space 
Telescope spectroscopy (D. Riechers et al., in prep.). 
Thus, IR pumping is likely substantially less effective 
for CN (relative to HCN, HCO+ and HNC), which may 
contrast^ the relatively high CN( 7V=4— >»3) luminosity 
suggested by Guelin et al. (2007), depending on the 
actual CN abundance. However, luminous CN(7V=3^2) 
emission was detected toward another distant lensed 
galax y, the Cloverleaf quasar (z=2.56; Riechers et al. 
l2QQ7bh . 

If CN and/or HC3N transitions close to those blended 
with HNC(J=5^4) and HCO+(J=5^4) are found to 
be much fainter than expected based on the V estimates 
by Guelin et al. (CN) and above (HC3N), a more com- 
plex scenario would be required to explain the different 
HCN, HCO+, and HNC line ratios in the J=5^4 and 
6^5 transitions. In this case, it will become necessary 
to investigate the different efficiency with which the IR 
radiation field excites their bending modes in more de- 
tail, which may lead to high but different excitation for 
the rotational ladders of these molecules. 

5.2. HNC/HCN Ratio 

Recent studies of HCN and HNC emission in nearby 
infrared-luminous galaxies have revealed sources with 
high HNC/HCN ratios and high HNC excitation. Two 
scenarios were brought forward to explain these high ra- 
tios: IR pumping and/or increased abundances of HNC 
in the presence of X-ray dominated regions (XDRs, often 
found in regions impacted by emission from active galac- 
tic nuclei; Aalto et al. 2007b). On the one hand, we find 
that HCN(J=6^5)/HCO+(J=6^5) > 1, which is the 
oppos ite to w hat is expected in the XDR scenario (Aalto 
et al. 2007b|). However, we note that our LVG mod- 
els show that the HCN(J=6^5) line (and thus, likely 
also the HNC J=6^5 line) is only moderately optically 
thick, abundance effects thus cannot be ruled out. On 
the other hand, the HNC(J=6^5)/HCN(J=6^5) ra- 
tio of 0.50+0.15 in APM 08279+5255 also fits well into 
the IR pumping scenario: in nearby infrared-luminous 
galaxies, the gas is warm enough to efficiently pump 
HNC (through the 21.6 /im bending mode with an en- 
ergy level of hv/k=669K and Einstein A coefficient of 
A m =5.2s- 1 ), but not HCN (through the 14.0 /im bend- 
ing m ode hu/k=1027K and Atr=1.7s _1 ; Aalto et al. 
\mmh . leading to a high HNC/HCN V ratio in high- 
J transitions. In APM 08279+5255, the dust and gas 
is warm enough to also efficiently pump HCN at high 
rates, so the line ratio depends mostly on the fraction of 
the HNC or HCN-emitting gas that is exposed to the IR 
radiation field and optical depth effects (given that the 
fraction of L' in the J=6^5 transitions due to collisional 
excitation is likely small). This allows for a broad range 
of HNC/HCN ratios, and thus is consistent with the com- 
paratively low J=6^5 ratio, especially in combination 
with the high excitation of both molecular line spectral 
energy distributions (SLEDs). Similar arguments can be 

9 The high CO excitation in APM 082 79+52 55 can be explained 
by collisional excitation alone (WeiB et al. 2007); however, due to its 
higher critical density, collisional excitation of CN is substantially 
less effective than for CO. 
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made for the relative strength of HC0 + , which has its 
fundamental 12.1 /im bending mode at slightly higher en- 
ergy than HCN. 

5.3. Detectability of Mir-IR Pumping Lines 

Given the SED and brightness of the high- J rotational 
lines of HCN, HNC, and HCO + , pumping by mid-IR ro- 
vibrational transitions appears plausible. As the abun- 
dance of these molecules is high enough to yield bright 
rotational lines in emission, the question arises whether 
or not the pumping transitions themselves may be de- 
tectable (which would place better constraints on the 
abundance of these molecules). The 14.0 /im HCN fea- 
ture was detected in absorption in nearby (U)LIRGs 
(Lahuis et al. 120071 ). Indeed, sources like Arp220 or 
NGC4418 show deep absorption features (10-30% of the 
continuum flux), indicating high HCN abundances, and 
evidence for pumping of the rotational HCN lines. The 
v 2 = 1 lines of HCO+, HCN and HNC are at rest- 
frame 12.1, 14.0, and 21.6 /im, i.e., observed-frame 59.4, 
68.8, and 106.1 /im, which is within the wavelength range 
of the Photodetector Array Camera and Spectrometer 
(PACS) on boar d th e Herschel Space Observatory (e.g., 
Poglitsch et al. l2010l ). APM 08279+5255 has continuum 
fluxes of 5 11+5 1 and 951+228 mJy at 60 and 100 /im (Ir- 
win et al. Il998f ). Thus, absorption features that have a 
depth of 10% of the continuum flux would be detectable 
within a few hours with Herschel. 



6. CONCLUSIONS 



HNC, 



Based on CARMA observations of HCN, 
and HCO + (J=6^5) toward the z=3.91 quasar 
APM 08279+5255, evidence is consolidating that the 
dense, star-forming molecular gas in this source is 
substantially enhanced in brightness by IR pumping (on 
top of gravitational magnification of the source) E3 The 
unusual, hundreds of parsec scale, warm gas and dust in 
this source appears to harbor a strong IR radiation field 



that efficiently pumps the high-J transitions of HCN, 
HNC, and HCO + at rates well beyond those achieved 
by collisional excitation alone. At least in the case of 
HCN, this leads to moderately optically thick emission, 
which appears as 'super-thermal' between J=5^4 and 
6^5 (due to a combination of high excitation and rising 
optical depth). These findings are consistent with the 
picture that the bulk of the gas and dust in this source 
is situated in a compact, nuclear starburst, where both 
the highly active galactic nucleus and star formation 
contribute to the heating. 

Even though APM 08279+5255 is known to be an ex- 
treme system, radiative excitation of dense molecular gas 
tracers may also play a significant role in other high red- 
shift galaxies, like it does in nea rby infrared-luminous 
galaxies (e.g., Aalto et al. 2007b). This issue becomes 
increasingly important for the interpretation of observa- 
tions of higher- J (J > 2) transitions of dense gas tracers 
like HCN, HCO+ and HNC, which are redshifted into 
the millimeter observing windows at high z, and thus 
will be primarily targeted toward galaxies in the early 
universe by future facilities such as the Atacama Large 
Millimeter /submillimeter Array (ALMA). 
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APPENDIX 
FLUX CALIBRATION 

The observations presented here were taken over the course of four (HCN & HCO + )/six (HNC) months. It thus is 
necessary to correct the flux scale of individual tracks for the intrinsic variability of the phase calibrator. Given the 
relatively large number of tracks (14/23), this leads to a substantially improved flux calibration relative to the sparsely 
sampled, internal look-up table of MIRIAD. 

Figure [5] shows the variability of J0818+423, the phase calibrator, bootstrapped relative to Mars, Uranus and 3C84 
(top panel), measured over the first band in the LSB at ~104.6/105.3GHz. For illustration, a fit to the data is shown. 
Fluxes for tracks that were observed without a primary calibrator were estimated based on a weighted average of those 
of tracks that are adjacent in time. As 3C84 is a millimeter- variable source itself J3 the bottom panel shows its variability 
over the same time range in which this project was observed. Within the calibration errors, the phase calibrator shows 
a clear excess variability over that of 3C84 in the 2009 measurements (data points at (JD-2454500)>332). Also, 
the scatter of values in the 2009 measurements bootstrapped relative to the Mars planetary model is substantially 
larger than that in 2008 (typically observed under comparable weather conditions). Together with the fact that the 
fluxes bootstrapped relative to 3C84 agree well with those bootstrapped relative to planet models, this suggests that 
J0818+423 shows significant intrinsic flux variations over the observed time range. Thus, we obtain the best calibration 
by using and interpolating the individually bootstrapped fluxes for each track, rather than using an averaged value 

10 IR pumping may also explain the unusual excitation of low- 
ionization rest-frame ultraviolet absorption lines in proximate ab- 
sorbers along selected sight lines toward the active galactic nucleus 
(which may arise from diffuse gas in the galaxy; Srianand & Petit- 
jean [2003). 

11 See, e.g., SMA data archive; 

http : / / smal . sma . hawaii . edu/ callist/ . 
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Fig. 5. — Millimeter variability of the phase calibrator J0818+423. Top: Measured ~105GHz flux for the tracks in this project where 
a primary calibrator was observed, bootstrapped relative to the flux calibrators indicated in the bottom right corner. Error bars indicate 
the variance of individual measurements within a track. The black line indicates a fit to the data. Bottom: Measured 95-112 GHz fluxes 
of 3C84 over the same time range (CARMA archive; S. Schnee, priv. comm.), bootstrapped relative to the flux calibrators indicated in 
the bottom right corner. Error bars include the absolute calibration uncertainty. This comparison shows that the variability among tracks 
bootstrapped relative to 3C84 is not due to intrinsic short-term variability of 3C84 itself. 

for J0818+423. We conservatively estimate that this flux calibration leads to an overall accuracy of ~15%. 
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